Gravitational wave energy spectrum of hyperbolic encounters 
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The emission of gravitational waves is studied for a system of massive objects interacting on 
hyperbolic orbits within the quadrupole approximation following the work of Capozziello et al. 
Here we focus on the derivation of an analytic formula for the energy spectrum of the emitted 
waves. We checked numerically that our formula is in agreement with the two limiting cases for 
which results were already available: for the eccentricity e — 1, the parabolic case whose spectrum 
was computed by Berry and Gair and the large e limit with the formula given by Turner 0]. 
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I. INTRODUCTION 



Einstein predicted already in 1916 that accelerated 
masses should emit gravitational waves. The detection 
of such kind of waves would open a new window in the 
exploration of our universe. In the last years technology 
has improved very rapidly, and it is now believed 
that the precision we reached should enable the direct 
detection of gravitational waves in few years, both with 
ground based and space based detectors such as e.g. the 
proposed eLISA mission. It is, therefore, interesting to 
study the dynamics of typical systems and their emission 
of gravitational waves and in particular their frequency 
spectrum, in order to know at which wave-length range 
we should expect gravitational radiation. 



For the cases of binary systems or spinning black holes 
on circular and elliptical orbits the resulting energy spec- 
tra have already been well studied [1, [f|. The energy 
spectrum for parabolic encounters has been computed 
either by direct integration along unbound orbits Q or 
more recently by taking the limit of the Peters and Math- 
ews energy spectrum for eccentric Keplerian binaries Q ■ 

The emission of gravitational waves from a system of 
massive objects interacting on hyperbolic trajectories us- 
ing the quadrupole approximation has been studied by 
Capozziello et al. [l[ and analytic expressions for the 
total energy output derived. However, the energy spec- 
trum has been computed only for the large eccentricity 
(j>1) limit Q. In this paper we derive the energy spec- 
trum for hyperbolic encounters for all values e > 1 and 
we give an analytic expression for it in terms of Hankcl 
functions. We checked numerically that our result in the 
limit of e = 1 is in agreement with the one for parabolic 
encounters @ and for large eccentricities with the result 
given in Q. 



II. THEORETICAL FRAMEWORK 

Gravitational waves (GWs) arc solutions of the lin- 
earized field equations of General Relativity and the 
radiated power to leading order is given by Einstein's 
quadrupole formula, as follows 



P = 



G 
45c 5 



{Dij Da) 



(2.1) 



where we used as definition for the second moment ten- 
sors Mij := J T 00 XiXj d 3 x, and for the quadrupole 
3Mii — SuMt-k- Here and in the 



moment tensor D 



following dots denote time derivatives 1 . 

The quantity My depends on the trajectories of the in- 
volved masses, and can easily be computed for all type of 
Keplerian trajectories. To compute the power spectrum, 
i.e. the amplitude of radiated power per unit frequency, 
requires a Fourier transform of equation (|2.1[) , which is 
rather involved (for the elliptical case see e.g. Q), and 
we will derive it below for hyperbolic encounters. 

In Fig. 1 the geometry of an hyperbolic encounter is 
represented with the most important quantities we will 
use. Since we will compare our results with those of Q 
and [l| , it is important to note that not all these quanti- 
ties are independent from each other, and we will need to 
know the relations between them. Notice that we assume 
that the gravitational energy loss during the encounter 
is negligible and thus that the Keplerian hyperbolic tra- 
jectory is a good approximation of the orbit. Clearly, 
this assumption does depend on the mass ratio and on 
the distance of closest approach. In the considered cases, 
where the masses are similar, this holds very well. 
The eccentricity e of the hyperbola is (see e.g. |f|) 



e:=A 1 
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1 Note that often in the literature (e.g. [(|) we also find the no- 
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tation Qij : = 
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equation 112.11 1 reads then: 
( Qij Qij Y- Here and in the following we will use 
the notation given by [J and [jj. 
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where E = 7j/i«o (E is a conserved quantity for which 
we can take the energy at t = -co), Vq being the 
velocity of the incoming mass mi at infinity, the angular 
momentum L = fxbvo, the impact parameter b, the 



reduced mass /i 



, the total mass m := m\+m2, 



TOi +7719 

and the parameter a := Gm^i 



We notice that the orbit is characterized by only four 
quantities given as initial conditions: vq, b and mi. 2. All 
the other parameters can be expressed as functions of 
this fundamental set. We can for instance rewrite the 
eccentricity as e = \JT 
the semi-major axis a = 



- vfjb 2 I G 2 m 2 = e(vo,b, m), or 
a/[j,v 2 , the angle at periastron 
through cos^o = — l/ e > or the radius at periastron as 



Gm 



(e — 1), and so on. 




present these computations, whose results we will then 
use. 

First we compute the energy and angular momentum 
of a body within a gravitational potential $(r). In the 
plane of the orbit the velocity can be written in terms of 
a tangent and a perpendicular component 



v r r 



tp, where v r = — 



dr 
di' 



~dt 



where vectors are represented by bold symbols. Thus the 
total energy per unit mass of the system and the angular 
momentum can be written as 



T 2 dtp 

L := r x v =r —— . 



dt 



Putting these equations together, using the substitution 
u := 1/r with r 2 = L/p> and rearranging, we get 



2E u 2 n 2$ f du\ , 2$ 

- ~ 2+u + ^ = U +u + i? 



L 2 tp 2 



FIG. 1: The geometry of an hyperbolic encounter and its most 
important parameters. 



Setting the angle of the incident body to <p = at 
initial time t = —00, the radius of the trajectory as a 
function of the angle and as a function of time is given 
by 



r(p) 



a(e 2 -l) 



1 



r(0 



cosh £ — 1) 



e cos{<p - ipo) 

with the time parametrized by £ through the relation 
= \p^f (e sinh£ — £), where £ goes from —00 to 

+00. 

Expressing this in Cartesian coordinates in the orbital 
plane, we finally get the equations for hyperbolic trajec- 
tories 



x (0 = a(s — cosh£) , 
y(Q — a Ve 2 — 1 sinh£ 



(2.3a) 
(2.3b) 



III. POWER SPECTRUM OF GRAVITATIONAL 
WAVES FROM HYPERBOLIC PATHS 

A. Power emitted per unit angle 

In [l[ the computation of the power emitted as a func- 
tion of the angle, as well as the total energy emitted by 
the system has been already carried out. Here, we briefly 



Since E and L are conserved quantities, the derivative of 
the last expression with respect to u gives: 

d 2 u 1 d$ d 2 u Gm 2 

This is an inhomogeneous linear differential equation of 
second order, which has the following solution 



u(<p) = B cos(</? — tpo) + 



Gm 2 



and substituting back to r we have r = B L sin(<p — (^0). 
B is a constant depending on the initial conditions and 
ipo is the polar angle corresponding to the periastron dis- 
tance, i.e. the distance of closest approach between the 
two interacting bodies (see Fig. 1). Its relation to the 
eccentricity is given by e = — 1 / cos po . 

Using the initial condition for the velocity, and the 
standard procedure to reduce the two-body problem to 
a single reduced mass particle moving in a gravitational 
field generated by the total mass, we see that the orbit 
of this reduced mass particle reads 



~(<P) = 



b sin po 



cos(p> — pa) — COS(/?o 



(3.1) 



In order to compute the emitted power P, given by 
the quadrupolc formula, as a function of the angle, it 
is convenient to rewrite the Cartesian coordinates Xi in 
spherical coordinates (the plane of the orbit corresponds 



3 



to -d 



x — r cos ip sin •§ = r cos if , 
y = r sin if sin 1} = r sin ip , 



3 = r cos i9 = . 
For the second momenta tensors My we get accordingly 



A/11 = /-i£ 2 = fJ-r 2 cos 2 y> , 

A/22 = MJ/ 2 = A"" 2 sm2 f ) 

-A/12 = M^'y = A 47 " 2 cos </? sm V 3 = -^21 i 

A/32 = A/23 = A/33 = A/13 = A/31 = 



(3.2) 



The term DijDij = J2i j DijDij in the expression for the 
radiated power can now be simplified to give 



<DijDij> = 6 <M n 



• ■•2 

A/22 



3A/ 



12 



A/llA/22> 



In order to compute this value explicitly, keeping if as a 
variable instead of t, we have to transform derivatives in 
time in derivatives in tp and r. This yields 



P(if) 



6, ,2 



32G£> 

45 c 5 



7(v.Vo) , 



(3.3) 



with 



sm((p - f ) 4 sin(|) 4 



150 + 72 cos(2(^ ) 



tan(< / 9 ) 2 sin(( / 9 ) 6 
66cos(2(</?o — ip)) — 144 (cos(2<^ — f) — cos(ip)) 



which is the result obtained in [l|. In Fig. 2 we plot the 
radiated power P as a function of the angle. 




FIG. 2: Radiated power as a function of the angle during an 
hyperbolic encounter, for e = 3, i.e. ipo — 0.6tt, according to 
the relation e = — 1/ cosi^o- 



Note that P depends on Dij, which also depends on the 2; 
present in the integral definition of Afy . 



Eq. (|3.3[) can also be written as follows, using the 
Schwarzschild radius r s := 2Gto/c 2 , and L = vq b 



P = 



dE 



' 45 c 3 b 2 



f(<P,<Po) 



Setting t = at pcriastron, the total energy radiated in 
GWs by the system during the interaction is given by 



AE 



00 


dE 


I 


dt 



dt . 



(3.4) 



Since we know ^ as function of if rather than t, we 
perform a variable change in the integration and get 



AE 



^r s vln r V2 sin 2 if q f (if, if Q ) 
45 c 3 b J [cos(f - f ) ~ cos f] 2 



dif , 



which can be evaluated taking if\ = as initial angle, 
and if 2 = 'Zfo as the final angle 



iJ5 = 32G^| 
bc° 



(3.5) 



with 
F(fo) 



720 tan 2 ipo sin 4 ip 



x [2628v?o + 2328^0 cos 2if 



+144(^0 cosAifo - 1948 81112(^0 - 301sin4^ 



This means that the total radiated energy of the sys- 
tem can be determined knowing the parameters b and Vq, 
and of course the mass \i. 

In deriving the above equation the implicit assumption 
has been made, that the energy loss doesn't change the 
path of the body in the gravitational field, which is rea- 
sonable since the emitted energy in form of gravitational 
waves is rather small. As stated above, this holds very 
well for the cases we consider here, where the mass ratio 
is about 1. 



B. Power spectrum 

We compute now P(uj), the Fourier transform of P(i), 
which describes the distribution of the amplitude of the 
power emitted in form of gravitational waves depending 
on the frequency. In the following we use the convention 



fit) e 



dt 



In [2| and [8j some hints are given when solving the anal- 
ogous problem in electrodynamics. The crucial idea is 
to use Parseval's theorem on the integration of Fourier 
transforms first, and then to express some quantities in 
terms of Hankel functions. This allows to compute in an 
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easier way the Fourier transform of P(t), for which we 
use the expression given in eq. (|2.ip 



AE = J P(t)dt = J P(lu)<:\lj = 7 / < Dij(t) Dij(t) >dt = 
7 / {\Du{uj)\ 2 + \D 2 2{io)\ 2 + 2\D 12 {uj)\ 2 + |S 33 (w)| 2 ) duj , 



where Dij{u) represents the Fourier transform of Dij(t), 
and we introduced the constant 7 := — -^p- 

It is easy to see that the last equation represents the to- 
tal amount of energy dissipated in the encounter. There- 
fore, the integrand in the last line has to be equal to the 
power dissipated per unit frequency P{oj), i.e. : 

p(u) = 7 (|i)iiHI 2 +|ii22HI 2 +2|i)i 2 (w)| 2 +|ij33MI 2 ) • 

_ (3.6) 

As next, we need to compute the Dij(uj), to square 
their norm and sum them together in order to get the 
power spectrum. Note that first we compute the Fourier 
transform of Dij(t), then we differentiate three times, 
and only at the end we take the square of their norm. 
Computing the Dij explicitly - keeping in mind that we 
use the time parametrization = \i a 3 jot (e sinh£ — 
£) - we get: 



D n (t) 

D 22 (t) 

D 12 (t) 



a 2 m 



2 

3 m a 2 



((3-e 2 )cosh2£-8ecosh£) , 
(4 e cosh £ + (2 e 2 - 3) cosh 2£) 
— (4 e cosh £ + e 2 cosh 2£) , 
\/e 2 - 1 (2 £ sinh £ - sinh 2£) 



(3.7) 



The Fourier transform of the third derivatives of Dij (t) 
is given by 



(3.8) 



thus we have just to compute Dij(w). To compute the 
Fourier transforms we can closely follow the calculations 
performed e.g. in where the similar problem in elec- 
trodynamics of the emitted power spectrum for scatter- 
ing charged particles on hyperbolic orbits is treated. In 
particular the following Fourier transforms are used 



sinh£ 



-Hg\ive) , cosh£ = -*H$\ive) , (3.9) 



cos 



where \x) is the Hankel function of the first kind of 
order a, defined as Ja(x) + iYa(x), with Ja(x),Ya(x) the 
Bessel functions of first and second kind, respectively; 



(3.10) 



and where the frequency v is defined as v := lu 
(the parameters /z, a and a as defined previously). 
Taking the above results for (t) we get 3 



D u (lo) 
-D22M 
-D33M 
D^ 2 {to) 



[16£HM'(iV£) + (£* -3)H£"(iV£/2)} , 



4w 
a rw, 

4w 
a 2 mi: 

4uu ^' 
3 a 2 mir 

4ue 



[(3-2£ 2 )H^'(w£/2)-8£H^'(iu£)} 

[8£H^'(W£)+£ 2 H^'(iP£/2)\ , 

y/^T [H$ (ii/e/2) - 4 £ H${ive)] 



Inserting this result into eq. (|3.6p , using eq. (|3.8[) , we get 
the power spectrum of the gravitational wave emission for 
hyperbolic encounters 



P(w) = - 



720 c 5 



where the function F £ (ui) is 



u 4 F e (u) , 



(3.11) 



F e (u) =\[lQ£H£»{iv£) + (e 2 3) H^(iu£/2)}\ 2 + 
I [(3 - 2£ 2 )H£ )f (w£/2) -8£H£ ) '(iv£)}\ 2 + 
\[8eH^'^s)+e 2 H^'^s/2)]\ 2 + 
9(e 2 - 1) 



\[H^{w£j2)-4,£H^{iv£)r . 

c 

In Fig. 3 the function uj 4 F e (u>) is plotted for some 
some values of e: this is the frequency power spectrum 
of gravitational radiation emitted by an hyperbolic en- 
counter. Unfortunately the expression for F e (uj) is rather 
complicated and we could not find an analytical way to 
simplify it. We thus made some numerical tests to check 
its validity and clearly the integral of (|3.11[) has to be 
equal to AE in (|3.5[) . which was obtained by integrating 
over the power emitted per unit frequency, i.e.: 

/ P(uj)duj = AE. (3.12) 
Jo 

We have checked the validity of this equality for dif- 
ferent sets of values, comparable to those used in []), e.g. 
b = 1AU, wo = 200 km/s, and mi, 2 = 1-4 M , or similar. 



3 Notice that we obtain also constant terms in these four equations. 
However, they can be dropped, since we can freely change the 
origin of coordinates while keeping invariant the description of 
the quadrupole radiation, as explained e.g. in @| in section 3.3.5. 
Note also that in fact these terms would in any case vanish for 
merely mathematical reasons, since the Fourier transform would 
multiply them with a Dirac delta function S(uj) and a factor 
(— iuj) :i from the third derivative. This expression vanishes for 
all uj ^ because of the 5(a)), and for uj = because of the 
multiplying factor (— iu)) 3 . 
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FIG. 3: The frequency power spectrum of gravitational radi- 
ation emitted by an hyperbolic encounter. On the x-axis we 
have the angular frequency ui expressed in mHz units, whereas 
on the t/-axis the amplitude of P{ui) is normalized to the max- 
imum value of the e ~ 2.5 case. These are the expected emis- 
sions generated by a system of two supermassive black holes 
with m = 10 7 A/q, impact parameter b = 10 AU, and differ- 
ent relative velocities. With lower velocities the interactions 
are stronger and the eccentricity decreases. These spectra, in 
order from the highest to the lowest, represent systems with 
wo = 3.4 x 10 7 m/s (e ~ 2.5), v = 3.5 x 10 7 m/s (e ~ 3), 
vo = 3.6 x 10 7 m/s (s ~ 3.1), v = 3.75 x 10 7 m/s (e ~ 3.4), 
v = 4 x 10 7 m/s (e ~ 3.8), v = 4.5 x 10 7 m/s (e ~ 4.7), 
respectively. In particular the case with e ~ 3 (plotted with 
the dashed line) is discussed in the conclusions. As one can 
see, for higher eccentricities the peak frequency slowly de- 
creases. This is only true for values of vo up to ~ 6 x 10 7 
m/s, whereas above it increases again. Moreover, decreasing 
the mass or increasing the impact parameter changes the ec- 
centricity as well. We should be able to detect incoming waves 
in that range e.g. with eLISA, since the peak at ~ 0.2 mHz 
fits in its observable band. For a more detailed discussion see 
Sec. IV. and e.g. [§]. 



For all of these sets we got agreement within numerical 
accuracy. 

More interesting is the case where the eccentricity ap- 
proaches e = 1. According to cq. (|2.2p this is the case 
e.g. with the set of initial conditions b — 2 AU, vq = 6.4 
km/s and mip = 1.4 M©. Since this is a limit case for a 
parabolic trajectory, we can directly compare our result 
with the one studied by 0, and indeed they coincide, 
within numerical accuracy. For a discussion about the 
feasibility of an analytical comparison sec Appendix B. 



interaction is written in Q as: 

G 7/2 m l/2 m 2 m 2 



AE 



15 C 5 7/2 



9{e) 



(3.13) 



where r m i n is the radius at pcriastron, and the enhance- 
ment factor g(e) turns out to be: 

24 arccos (^) (l + |f e 2 + f e 4 ) + V^~T{ ^ + gge 2 ) 

(e + iy/ 2 

Clearly this expression is equivalent to our result (|3.5[) 
for all values of e > 1. Indeed, using the following rela- 
tions for the radius r m i n and the angle tpo at periastron: 



-1 



COS (fo 



sin tpo 



1 — COS ifio v { 



both expressions can be written as functions of only three 
parameters describing the encounter: e, m and vq, and 
it is then straightforward to see that 

AE Q = AE T , 

where AEq denotes the energy computed according to 
eq. p.5p and AEt the one according to eq. (|3.13[) . 

Following Q again, we find out that in the limit for 
large e, the g(e) factor can be simplified and written as: 



(3.14) 



which also agrees with the result of Wagoner and Will 
[loj . This yields a simple form for the energy emitted 
during the path: 



AE T = —— 5(e) , 



15 c 5 r 7 / 2 

min 



(3.15) 



That leads then to the formula for the energy spectrum 
valid for large e: 



P(a) = 



G 7/2 m l/2 TO 2 TO 2 g 

7 / 2 15tt 



c° r 



x |l2[(T 2 if 2 (a) 



mm 

-aKi{a)] 2 +3[2a 2 K 1 (a) + <jK {a)} 2 +o- 2 A 2 (cr)| , 

(3.16) 

where Ka(x) are the modified Besscl functions of the 
second kind, a is the frequency rescaled in terms of 
the characteristic time scale of the gravitational wave r 
[=(periastron distance)/(periastron velocity)], a = lot. 



C. The limit for e > 1 

As next we turn to the large e limit and compare our 
result with the one given in 0] and [Toj]. The expres- 
sion for the total energy emitted during an hyperbolic 



Comparing our total energy from the quadrupole ap- 
proximation, eq. (|3.5[) . with the expression for the energy 
AE T (|3~T5]) by 0] with the simplified factor ([3714]) valid 
in the large e limit, we see that they coincide for large 
eccentricities, having e.g. a 1% difference after e = 100, 
and a 5% difference after e = 20 as shown in Fig. 4. 
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In fact we see that the behavior of the variation goes 
as (AEq — AE T ) I AE T oc 1/e, confirming the fact that 
for the parabolic limit AEq = 2 AEt, so that one would 
underestimate the energy emitted by a factor of 2 taking 
this approximation Q. 



o.io - 

h 0.08 - 
<1 

«o 0.06 - 
<1 

< 0.04 - 

0.02 - 

0.00 t 




50 100 150 



FIG. 4: The converging behavior of the total energy emit- 
ted in an hyperbolic encounter according to the quadrupole 
approximation AEq, eq. (|3.5[) , towards the result of 0], eq. 
IpTTS)) . The plotted line is (A£q - AE T )/AE T as function 
of the eccentricity e. 



IV. CONCLUSIONS 

Short gravitational wave burst-like signals are expected 
in the data stream of detectors. Although these signals 
will likely be too short to allow us to measure the pa- 
rameters of the emitting system accurately, the results 
presented in this paper could be used to get a rough es- 
timate of these parameters, by observing the position of 
the peak, the amount of energy released and the timescale 
of the interaction. 

Given the knowledge of the power spectrum we can 
easily see which kind of hyperbolic encounters could gen- 
erate gravitational waves detectable e.g. with eLISA, ad- 
vanced LIGO or advanced VIRGO. Measurements from 
unbound interactions with ground-based detectors could 
in principle be possible, though the energy emitted at e.g. 
±200 Hz is below the minimum threshold for advanced 
LIGO or advanced VIRGO, making detections unlikely 
but not impossible. The space-based interferometer in- 
stead is expected to cover frequencies ranging from 0.03 
mHz up to 1 Hz (see e.g. [9(), where the interactions 
could release more energy. 

An unbounded collision between two intermediate- 
mass black holes, let's say of 10 3 M Q each, with an en- 



counter velocity of 2000 km/s at a distance of 1 AU, 
would generate, according to our eq. p. lip , a frequency 
spectrum with peak around 0.04 mHz, with 80% of the 
emission in the range between 0.01 and 0.07 mHz, i.e. in 
the lower range limit of eLISA. Another possible exam- 
ple of measurable impact would be an encounter between 
two supcrmassivc black holes with mass, e.g., compara- 
ble to the expected mass of Sagittarius A*, the black hole 
believed to be at the center of our galaxy, i.e. ~ 10 7 M Q . 
With a distance of some AU, and a high velocity (we 
want to exclude the bounded case) of tens of thousands 
km/s, such a collision would generate an energy spectrum 
with peak at ~ 0.2 mHz with 80% between 0.03 and 0.37 
mHz, thus in the observable range of eLISA. (Its energy 
spectrum is plotted with a dashed line in Fig. 3.) 

Interestingly, the time window of such events is enough 
to allow measurements. Indeed, for encounters up to e ~ 
5 with peak in the mHz to Hz regime, we are in the time 
scale of 1 day, if we choose the cutoff of the interaction 
at an angle of ip = 3/4ifQ, i.e. where the path starts to 
approach significantly the asymptote (see Fig. 1). For 
instance, for the two examples above, the interactions 
would last about 54h and 9h, respectively. Estimates for 
the rate of such events have been considered e.g. in [l3| . 
They consider e.g. typical compact stellar cluster around 
the Galactic Center, and expect an event rate of 10~ 3 up 
to unity per year, depending on the radius of the object 
and the amount of such clusters in the near region. 

After a discussion with L. Blanchet we realized that 
there is a possibility to treat the problem in an alter- 
native way. Starting from the Keplerian equations of 
motion, one could take the solution of Peters and Math- 
ews JH, bring the argument onto the imaginary axis and 
- making use of equations (9.6.2-9.6.4) in Abramowitz & 
Stegun [11| - find the energy spectrum in terms of Hankcl 
functions. 

We believe that with the wave-form found here one 
should be able to classify the different encounters de- 
pending on the detected shape, and therefore get a better 
insight into the map of our galaxy or the near universe. 
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APPENDIX 
A. On the Fourier Transform of sinh£ and cosh£ 

In section 1111 Bl wc used the relation (|3.9[) without any proof. Since we didn't find any reference where the proof is 
shown explicitly, we will show in the following where these two relations arise from. In Landau and Lifshitz the 
equations we used in the cited section are: 

St = -—H£ ) (ive) , c^sh| = --H$'(ive) (A.1-A.2) 

but in fact the general relations we want to show - [3], §70, equation (70.15) - are: 



m = a -^^H2){ive) , m = -H^'^e) (A.3-A.4) 
I. Fourier Transform of y(t). We know that the time can also be parametrized with £ using the transformation: 



= \l — — (e sinh£ — £) which allows to write y{£) = a \/ e 2 — 1 sinh£ 
We can now start with the computation of the Fourier transform: 



J-oo J-oo <K J -00 M <K 

IU J-oo dt, J~oo 

aV^r f°° , , , 1 1 



/oo 1 1 

(cosh£-- + -)e- iw * ( « ) d£ = 

/ (e cosh£-l)e- iwt ^d£+ / 

J —OQ J — 



e -iut{&) d £ 



We now compute separately A and 6: 



.^d^Vfia 3 V fxa 3 J_ OCl iu d£ 



Hence we have: 



2/M = 



O + iTrf&^ii/e) 



aVe 2 - Itt (i),. . 



(A.5) 



□ 

where: (* : ) follows from partial integration and the boundaries vanish; we have the same situation for (* 2 ); and in 
(* 3 ) we used the defining relation of the Hankel functions. 
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II. Fourier Transform of x(t). We now show equation (A. 4) in a similar way: again with the time parametrization 
through £ we can write: 



x (0 = a(£ — cosh£) 



Its Fourier transform can be taken as follows: 



x(u) : = / x(t) e - *"* dt = / x(£) e 



,t(5) 



dm 



1 d 



d£= / x(£) — — t e - iwt( -® d£ 



/OO J /*00 

(e-cosh0 4:e- iut «^= — / sinh£ 
-oo ^ 



—a 
ja; 

-a 1 
ito 2 
-a 1 
2 

-a 1 
2 



dC 



OO 



ive sinh 



g — li/£sinh£+(ii/+l)£ ^ _ 



-ii/e sinh 1)£ 



a7r 1 
~uT 2 



Hence, using the formula for the derivative given in (|3.10p we get: 

x (u) = — H\J (ive) 

UJ 



(A.6) 
□ 



B. Analytical comparison of the parabolic limit 

In section Till Bl we discuss the parabolic limit of our result for the energy spectrum, where e — > 1. A numerical 
check of its validity with respect to the previous known quantity, in 0] , is quite straightforward. Nevertheless it would 
be interesting to find an analytical agreement between these two formulae. 

Taking that limit the last part of our equations drops away, and only three terms remain. These can be greatly 
simplified, since in the e = 1 limit the second and third arc similar and F £ (ui) becomes: 

F 8 (w) = [16 H^' (ive) -2H^'(ive/2)} 2 + 128 [H$'(ive) ] 2 + 2 [H^'(ius/2) ] 2 (B.l) 

Following the book by Abramowitz & Stegun [ll| and the work done in Q we see two facts. First, the frequency in 
the order of these functions goes as (1 — e)~ 3 / 2 (see eq. (10) in Q), and therefore in the limit e — > 1, v goes to infinity. 
Second, we note that the order is also in the argument, therefore when v — > oo we can write the Hankel functions in 
terms of Airy functions Ai(x), e.g. eq. (9.3.45) in (Tlj : 

afiV.) ~ coe* x I^HO g ^ + A^ffO^ (B , 2) 

fc=0 fc=0 J 

where f is a function of e defined as §C 3 ^ 2 = l n 1+v,1 ~ ea — y/l — e 2 , see eq. (9.3.38) in If we only take the first 
order of the series, we have for the derivatives a term in Ai(x) and one in Ai(x). Again using fllj | we find eq. (10.4.26 
- 10.4.31) which express these terms in form of modified Bessel functions of the second kind: 

tf±i/ 3 (0 = irV3AAi(*) , #±2/3(0 = -*V3/zAi'(z) (B.3) 
where z = (§C) 2 ^ 3 - We see now that equation (B.l) has the very same structure as eq. (24) in Q: 

1(f) = [ 8B(f) 2.4(/) ] 2 + C [A(f) ] 2 + D [A(f) ] 2 (B.4) 

where A(f) and B(f) are also in terms of modified Bcsscl functions of the second kind of order ±1/3, ±2/3. 

At this point with some cumbersome algebra one should find that also all the coefficients of (B.2) and (B.4) agree. 
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